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ABSTRACT 


A  summary  is  presented  of  the  derivation  of  the  equation  of  notion 
of  the  bubble  due  to  an  underwater  spherical  explosion.  Migration  and 
surface  effects  are  included.  The  formulae  for  the  first  maximum  radius 
and  the  first  oscillation  period,  including  surface  effects,  are  given. 

From  the  experimentally  determined  bubble  period  data  and  the 
theoretical  period  formulas 

Tg  -  KW1’/3Z“5/6|l.O  ♦  KW1/3Z“1AffAD«-B)j| 

the  energy  partitions  of  three  explosives,  Pentolite,  BOX  ♦  Alum  (50/£0), 
and  2H2  ♦  Og  are  calculated. 
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SUMMARY 


The  Problem 

To  compare  underwater  explosions  producing  steam  bubbles  and  non¬ 
condensable  gas  bubbles  on  the  basis  of  their  computed  energy  partition 
values.  To  determine  which  of  three  selected  chemical  explosives  best 
simulate  an  underwater  nuclear  explosion. 


Findings 

Using  Friedman's  formula  for  the  bubble  oscillation  period  and 
the  experimentally  determined  bubble  periods  from  Hydra  studies,  the 
energy  partition  values  (i.e.,  fraction  of  charge  energy  left  for 
bubble  oscillation  after  shock  passage)  of  three  explosives  were  found 
to  be: 


Explosive 
Fentolite 
REK  +  Alum 

2H2  +  02 


r 

0.45  +  0.12 
0.50  +  0.14 
0.41  +  0.11 


Since  the  minimum  error  variation  in  energy  partition  is  much 
greater  than  the  variation  due  to  type  of  explosive,  it  makes  little 
difference  which  explosive  (steam  or  non-condensable  gas)  is  selected 
to  simulate  the  nuclear  case. 
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INTRODUCTION 


Undervater  nuclear  explosions  form  steam  bubbles  whose  dynamics  de¬ 
termine  the  early  dispersion  of  the  radioactivity  produced.  To  simulate 
nuclear  explosions,  the  Hydra  program  had  been  using  chemical  high  ex¬ 
plosives  (e.g.  Pentolite)  which  form  non-condensable  gas  bubbles.  In 
an  attempt  to  more  accurately  simulate  these  detonations,  two  steam 
producing  explosives  were  developed:  an  equal  mixture  by  weight  of  HEK  + 
Alum  (i.e.,  AlNHj^SOj^g-lSlHgO),  and  a  pressurized  stoichiometric  mixture 
of  hydrogen  and  oxygen.  Using  these  two  explosives  and  Pentolite,  a 
series  of  experiments  was  undertaken  during  the  summer  of  i960  to  deter¬ 
mine  the  variation  of  bubble  parameters  between  chemical  explosives  pro¬ 
ducing  condensable  steam  bubbles  and  non-condensable  gas  bubbles. 

The  bubble  parameter  usually  compared  is  the  energy  partition  value 
(i.e.,  the  fraction  of  charge  energy  left  for  bubble  oscillation  after 
shock  passage).  This  is  best  determined  from  the  theoretical  formula 
for  the  bubble  oscillation  period  using  experimental  data. 

This  report  reviews  the  bubble  dynamics  theory  that  has  been  deve¬ 
loped  over  the  past  several  decades  including  surface  effects.  Using 
these  theories  and  the  bubble  period  data  from  Hydra  studies,  the  energy 
partition  values  for  selected  explosives  of  interest  have  been  determined. 

In  the  literature,  the  first  order  approximation,  relating  the 
bubble  period  constant  K  to  the  cube  root  of  the  bubble  energy  is  well 
known.  However,  the  second  order  effects  of  fixed  and  free  surfaces  on 
the  bubble  period  are  generally  ignored.  Where  these  effects  are 
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mentioned,  the  approaches  are  difficult  to  follow  because  of  inconsist¬ 
encies  or  incompleteness .  After  a  critical  review  of  the  literature, 

•  Friedman's  papers  were  found  to  be  the  first  and  only  rigorous  qualitative 
and  quantitative  determination  of  surface  effects. 

Because  no  papers  were  found  to  use  Friedman's  method  to  its  full 
advantage,  and  because  Friedman's  original  work  is  more  detailed  than 
necessary  for  the  field  engineer,  the  following  section  on  bubble 
dynamics  was  undertaken  with  effort  to  present  a  comprehensive  and  co¬ 
herent  outline  of  the  equations  of  motion  using  a  consistent  system  of 
symbols  in  a  concise  package  permitting  bubble  period  prediction. 


SECTION  1 


DYNAMICS 


1.1  EQUATION  OF  MOTION 

The  earliest  derivation  of  the  equation  of  motion  of  the  gas  bubble 
Is  that  of  ftamaauer.3-  He  used  the  energy  balance: 

E  »  E  +  WD  +  KE  (l.l) 

o 

which  translated  to  words  says  that  the  initial  Internal  bubble  energy 
Eq,  equals  the  internal  energy  at  a  later  time  E,  plus  the  wort  of  dis¬ 
placement  against  the  water  WD,  plus  the  Xlnetlc  energy  imparted  to  the 
surrounding  water  at  any  later  time  XE.  Here 

Eo  ■  PoV<>'-1> 

E  -  FV/(r-l) 

WD  -  Ph(V  -  VQ) 

KE  -  2*  f>  A3A2 

Assuming  the  bubble  expands  adlabatlcally,  then 

FT7  -  P0VQ7  (1.2) 

Note:  See  Appendix  7  for  definition  of  symbols. 
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and  Eq.  1*1  can  be  rewritten  as 


P0/(r-l)  [uo/A)3  -  (A^A)3*]  ♦  Ph  ^A)3  -  lj  -  3/A2/2  (1.3) 

In  the  existing  literature  these  equations  are  usually  written  in  dimen¬ 
sionless  variables  a  and  t,  by  using  the  unit  of  length  L,  and  the  unit 
of  tine  C: 

L  -  (3rgw/Uy>gZ)1/3  (l.U) 

C  -  (3/2gZ)l/2  L  (1.5) 

Thus  a  ■  A/L,  t  *  T/C,  and  a  "  CA/L.  Subscripts  0,  1,  and  2,  refer  to 
variables  at  detonation,  first  maximum,  and  first  minimum. 

Since  EQ  ■  rQW  -  L3UnPh/3,  then 

P0/(r-D  -  Pha0“3  -  Phk-V(r-1)  (1.6) 

where  k  is  a  convenient  notation  and  is  defined  by 

k  -  a03(Y-1)  -  (UiyqgAo3V3rC^)(Y“:L)  (1.7) 


and  Eq.  1.3  now  reads 


1.0  -  ka 


-3(Y-1) 


♦  &  -  a^ 


3*2 

a  a 


Since  P^Y-l)  »  Ph'  that  is,  1.0  »  aQ3,  the  aQ3  term  is  usually 
dropped  and  the  equation  is  written t 


1.0 


ka 


-3(Y-1) 


3*2 

arar 


(1.8) 


This  approach  does  not  include  the  effects  of  fixed  and  free  surfaces. 


k 


and  assumes  no  migration 


1.2  MIGRATION 

During  the  first  half  oscillation,  the  migration  is  found  from  ex¬ 
periments  to  be  very  small.  However,  including  migration,  the  kinetic 
energy  term  would  be2 

Kb  ■  2  Y>A3  (A2  +  32/6) 

Equation  l.B  would  now  read 

1.0  -  ka“3(Y"1J  -  a3  -  a3(a2  ♦  b2/6) 

2 

The  migration  is  usually  written  in  terms  of  a  non-dimensional  vertical 
momentum,  s.  Since 

s  -  a3S/3  (1.9) 

then  the  equation  of  motion,  including  migration,  is 

1.0  -  ka'3^"1^  -  a3  -  3s2/2a3  -  a3a2  (1.10) 

At  the  maximum  and  minimum  radius,  a  *  0.  Thus  a^  and  Sg  are  the 
largest  and  smallest  roots  of 

1.0  -  ka"3^'1^  -  a3  -  3s2/2a3  -  0 

(Notet  aQ  is  also  a  root  of  the  equation  retaining  the  a  3  term  since 
sQ  -  0.) 

The  maximum  limit  of  ai  approaches  1..0  as  k  (l.e.  aQ)  and  s  approach 
zero.  Thus  a  good  approximation  to  tn-  maximum  bubble  radius  Is  given  by 
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(1.11) 


^  -  1.0  -  k/3  -  (3r-2)k2/9  -  s12/2 

(Reference  2  gives  3 3-^  /2.  This  is  a  typographical  error.) 

In  determining  the  oscillation  period,  the  momentum  term  is  negli¬ 
gible  (as  shown  by  calculations  in  Appendix  2),  and  therefore  is  ignored 
in  all  further  discussion.  Thus,  the  time  to  first  maximum  radius  is 
found  from  Eq.  1.8  to  be: 

ti  •  a3/2  (1.0  -  a3  -  ka"3^"1^)"1^2  da  (1.12) 

1'iaking  use  of  Friedman's  1  functions  (see  Appendix  1),  t^  *  1^.  Assum¬ 
ing  t2  “  2t^,  i.e.  assuming  the  integrations  from  aQ  to  a^  and  from  a^ 
to  a2  are  identical,  then 

T2  ■  Ct,  -  2CI1  (1.13) 

1.3  SURFACE  EFFECTS 

The  bubble  motion  is  affected  by  surfaces,  i.e.  the  bottom,  targets, 
and  the  air-water  interface.  Free  surfaces  repel  while  fixed  surfaces 
attract  the  bubble.  The  oscillation  period  is  reduced  by  the  free  sur¬ 
face  and  increased  by  the  fixed  surface.  These  effects  do  not  cancel; 
the  free  surface  tends  to  be  stronger.  Thus  the  equation  of  motion  must 
include  surface  corrections.  The  only  rigorous  surface  correction  method 
available  is  that  of  Friedman  in  reference  2  and  summarized  below.  A 

3 

comparison  of  Friedman's  final  equation  with  those  of  Kennard  is  also 
included. 

1.3*1  Friedman's  equation 

In  determining  the  effect  of  plane  (or  spherical)  surfaces,  the 
above  approaoh  to  the  equation  of  motion  is  valid  with  one  exception. 

The  kinetic  energy  of  the  water  is  found  using  a  velocity  potential 
function  which  is  evaluated  using  the  method  of  images.  This  method 
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is  detailed  in  reference  2.  The  resulting  kinetic  energy  expression, 
ignoring  second  and  third  order  terms  and  migration,  is 


M  -  2iy»A3A2(1.0  +  AF/(D+B)) 

where 

F  -  2xf (x)  -  ln2  (1.1U) 


x  -  (d-b)/(d+D)  (1.15) 

£(x)  -  2  n  (-1.0)nf(2n+l)2  - 
n>*0  I 

For  tabulated  values  of  x  versus  F(x),  see  Appendix  3- 
The  non-dimensional  equation  of  motion  is  now 


r 


(1.16) 


1.0  -  ka"3^"1^  -  a3 


a3  a2 


(1.17) 


Since 

order 


*  |l.O  +  aF/(d+b) 

!1.0  +  aF/(d+b)j1//2  -  1.0  +  aF/2(d+b)  + 
erms,  the  time  to  first  maximum  radius 

ti  -  J1  a3/2  (1  -  a3  -  ka"3^’1))"^2(l  ♦  aF/2(d+b))da  (1.18) 


...,  then  ignoring  higher 
is 


Or  in  terms  of  the  I  functions, 


hm  h*  y/2ta*b)* 


And  the  oscillation  period  is 


T 


2 


21-jC 


1.0  +  IgLF^I-^D+B) 


(1.19) 


If  the  explosive  were  detonated  at  a  depth  such  that  F  -  0  (i.e.  at 
D  •  2B),  then  Tg  ■  2CI1#  Since  this  must  agree  with  the  well  known 
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formula 


T  *  KW^V5/6,  (l.20) 

then  the  period  constant  K  is  defined  fcy 

K  -  2I1(3/2B)1^2(3rQA*^8)^3  (l.2l) 

For  convenience,  a  constant  is  defined: 

H  =  Ig(2g  (1.22) 


Now  Eq.  1.19  reads 

T  =  KM1/,3Z"5//6 


1.0  +  KW 


1/3, 


f1/3HF/(I>+B)j 


(1.23) 


This  is  the  form  of  Friedman's  theoretical  period  formula  which  will  be 
used  later  to  compute  energy  partition  values  (see  eq.  2.3). 

Next  consider  the  first  maximum  bubble  radius: 

Anmx  =  L(l*°  ’  "  (3 7-2)k2/9)  (1.24) 


Prom  Eq.  1.7,  k  does  not  equal  zero,  since  AQ  cannot  be  zero.  However, 
in  experiments,  the  error  in  the  measured  maximum  radius  is  much  greater 
than  the  second  order  correction  due  to  k.  Thus  =  L  must  agree 
with  the  well  known  formula 

A^  =  JW  1/3z”1/3  (1.25) 

This  defines  the  explosive's  radius  constant  J  to  be 

J  -  terfiAvg)1^3  U-26) 
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The  commonly  used  K/J  ratio  is  seen  to  be 

K/J  -  2I1(3/2g)1^2  (1.27) 

All  the  I  functions  depend  on  k  and  ft  but  from  the  graphs  of  I(k,y)»  in 
reference  2,  1^  >  0.7405.  Thus* 

*/ J  *  0.320  (1.28) 


2 

From  the  sane  set  of  graphs,  IjVig  2  0.92,  and 

H  <  1.26  (1.29) 

1.3.2  Kennard's  equation 

Kennard's  reports  of  the  early  191*0 's  often  give  the  formula 
T  -  uHh-S'b  (1.0  t  0.2Ajmuc/R) 

where  R  is  the  distance  from  the  dharge  to  the  fixed  or  free  surface* 

The  ♦  sign  is  associated  with  the  attractive  force  of  the  fixed  (bottom) 
surface,  and  the  -  sign  with  the  repulsive  force  of  the  free  (air-water) 
surface.  Thus  the  period  would  be 

T  -  Bf^V5/6  (1.0  ♦  0.2^  (D-B)/DB)  (1.30) 

Usually  Eq.  1.30  overcorrects  the  bottom  effect,  and  also  the  surface 
effect  in  deep  shots. 

In  reference  3,  Kennard  says  1.30  is  an  approximate  eqiirloal  ex¬ 
pression  and  that  the  theoretical  formula,  using  Friedman's  method,  is 

i .  (1.0  ♦  O.ZPMA^dHBj/DB)  (1.31) 
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where  x  is  defined  by  j£q.  1.15  and  values  of  P(x)  are  calculated  from 
reference  2.  Comparing  Kennard's  graph  of  P(x)  with  Friedman's  tabu¬ 
lation  of  F(x),  we  find 

P(x)  -  (1-x2)F/2  (1.32) 


Since  *  Le^,  Bq.  1.31  can  be  written  as 


KW^V^l.O  ♦  O.Ua1LF/(D+B) 


(1.33) 


This  agrees  with  Friedman's  expression  (Eq.  1.19)  if  lg/2Ii  * 

Thus  Kennard's  eq.1.33  is  a  special  case  of  Friedman's  eq,  1.23. 

A  choice  should  be  made  between  Friedman's  theoretical  period  eq. 
1.23  Involving  the  complicated  F(x)  function,  and  the  much  simpler  eq. 
1.30  of  Kennard.  Friedman’s  eq.  1.23  Is  preferred  not  only  because  of 
its  theoretical  foundation  but  also  because  It  better  fits  the  experi¬ 
mented  data  (see  page  13) • 
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SECTION  2 


ENERGY  PARTITION 

2.1  ENERGY  YIELD  OP  EXPLOSIVES 

Detonation  of  an  explosive  releases  an  amount  of  energy  (®f.  An 
amount  (l-r)SJtf  is  carried  off  by  the  shock  wave,  leaving  the  bubble 
with  rQW  for  the  first  oscillation. 

Reported  values  of  Q  for  any  one  explosive  are  varied.  For  example, 
for  the  standard  explosive,  TOT,  they  range  from  around  500  gram  calories 
per  gram  weight  to  over  1000.  Keeping  in  mind  that  the  accuracy  of  Q 
may  be  improved,  TOble  2.1  was  used  in  reducing  the  experimental  data 
of  this  report. 

TABIE  2.1 
Energy  Yields 


Explosive 

Q* 

(cal/gm) 

Q 

(ft-lb/lb) 

TOT 

1000 

1.40  lof 

Pentolite 

1220 

1.71  10? 

RDX  +  Alum  (50/50) 

1335 

1.87  ID? 

2R2*°2 

3794 

5*31  30 

Q  values  are  usually  reported  in  calories  per  gram  weight.  These 
have  been  converted  to  the  fps  system  using  the  conversion  factors: 

453*6  grams  =  1.0  pound 

1.0  gram  calorie  =  3*087  foot-pounds 


♦Private  cctmuni cation:  W.  V.  Perkins,  USNHDL. 
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The  Q  value  for  2Hg  ♦  Og  is  based  on  the  change  In  enthalpy. 

2.2  EVALUATION  OF  K  AND  J 

The  fraction  of  charge  energy  left  in  the  bubble  for  the  first 
oscillation,  that  is  r ,  can  be  found  using  Eqs.  1.21  or  1.26.  Thus 


( Uyg/3Q)  K3  ( 2g/3  )3/2/8l13 

(2.1) 

(bygMJ3 

(2.2) 

where  I-^  is  a  function  of  k  and  y,  Q  is  the  total  charge  energy  per 
unit  weight,  and  K  and  J  are  constants  of  the  explosive  to  be  deter¬ 
mined  from  experimental  data. 

In  section  1.3,  the  general  form  of  the  equation  for  the  oscil¬ 
lation  period  was  found  to  be  fiq.  1.23.  Written  as 

y  -  K  +  HX2z,  (2.3) 

where 

y  -  TZ^W1/3 

and 

z  -  w^3z'1/3f/(d+b), 

this  is  the  equation  for  a  straight  line  with  slope  HK2  and  intercept 
K.  Using  the  method  of  least  squares  (Appendix  5)  to  fit  the  explosive’s 
experimental  data  (for  varying  charge  depth,  pond  depth,  and  perhaps 
charge  weight),  the  values  of  K  and  H  can  be  confuted  for  the  explosive. 

Using  a  least  squares  fit,  it  is  possible  to  minimize  the  error  in 
y  or  in  z.  From  the  experimental  data  of  Appendix  U,  the  first  order 
variable  y  is  at  least  a  factor  of  100  times  the  second  order  variable 
z.  Therefore  the  following  values  of  K  and  H  were  computed  based  on 
minimizing  the  error  in  y» 
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TABLE  2.2 

First  and  Second  Order  Period  Constants 


Explosive 

K 

K' 

H 

Pentollte 

4.59 

4.30 

1.23 

REX  +  Alum  (50/50) 

4.89 

4.44 

1.32 

2H20  +  02 

6.50 

4.17 

1.50 

This  needs  seme  explanation.  The  K  'values  are  based  on  actual  charge 
weight  (except  for  HOC  +  Alum  where  0.8  pounds  was  used;  that  is,  the 
Alum  weight  was  excluded) .  However  If ,  Instead  of  the  actual  weight, 
the  equivalent  weight  TNT  (i.e.  QW/Q^)  is  used,  the  K*  values  com¬ 
pare  with  those  found  In  the  literature  which  vary  from  4.19  to  4.37 
for  TNT  and  Pentollte. 

From  Eq.  1.29,  the  largest  possible  value  of  H  is  1.26.  Thus 
except  for  Pentollte  the  above  H  values  are  too  high,  but  these  high 
values  can  be  explained  by  experimental  error  (see  Appendix  6).  If 
the  error  in  z  were  minimized,  the  H  values  would  be  still  higher. 

J  can  be  determined  from  Eq.  1.27  If  the  value  of  I,  (k,7)  is 
known.  The  ratio  of  specific  heats  of  the  bubble's  gas  is  not  always 
known,  but  it  is  standard  procedure  to  let 

7  ■  1.25  (2.4) 

Thus  to  determine  J,  k  must  be  evaluated.  But  freu  Eq.  1.7,  k  is  a 
function  of  the  energy  partition  r. 

2.3  ENEHX  PARTITION  VALUES 

As  shown  in  section  2.2,  the  energy  partition  can  be  determined 
from  Eq.  2.1  by  a  trial  and  error  method  since  r  ■  f(lj)  -  f(k)  ■  f(r). 
In  other  words,  it  is  a  circular  relationship.  To  simplify  things,  we 
solve  for  k(l^)  by  eliminating  r  between  Eqs.  1.7  and  1.21.  Thus 
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(2.5) 


k  -  ^(z/wKs/ag)3/2^^)3  j  (r  X) 

From  the  experiments,  the  charge  radii  in  feet  are:  Pentolite, 
0.1354;  REK  +  Alvun  (50/50),  0.1458;  2Hg  +  Og,  0.50.  Equation  2.5  gives 
one  relationship  between  k  and  1^.  Fran  reference  2,  the  graph  of  1^ 
versus  k  for  7  ■  1.25  gives  a  second  relationship.  The  intersection  of 
the  two  curves  gives  the  sought  for  values  of  k  and  1^.  J  is  found  from 
K  and  1^.  Knowing  k  yields  I^/lg,  and  thus  H  from  Eq.  1.22.  Finally, 
r  is  found  from  Eq.,  2.1.  The  results  are  below. 


TABLE  2.3.1 

Energy  Partition  Values  (7  -  5/4) 


Explosive 

k 

*1 

J 

J' 

H 

r 

Pentolite 

0.075  +  .001 

.7423 

14.32 

13.42 

1.23 

0.45 

RE0C  +  Alum 

0.081  +  .001 

.7420 

15.26 

13.86 

1.23 

0.50 

2Hg  +  Og 

0.213  +  .008 

.7410 

20.32 

33.03 

1.14 

o.4i 

where  the  +  values  are  due  to  the  depth  (Z)  variation. 

For  sea  water,  l.e.  pg  ■  64  lb/ft3,  the  respective  r  values  are 
0.46,  0.51  and  0.42. 

The  effect  of  increasing  7  from  5/4  to  4/3  is  shown  below. 

TABIE  2.3.2 

Energy  Partition  Values  (7  *  4/3) 


Explosive 

warn 

ri 

J 

J* 

H 

r 

Pentolite 

0.032  +  .001 

.7480 

14.21 

13.31 

1.24 

0.44 

REK  +  Alum 

0.035  +  .001 

.7482 

15.14 

13.74 

1.24 

0.49 

2Hg  +  02 

0.127  +  .006 

•7550 

19.94 

12.79 

1.17 

0.39 

Again,  for  sea  water,  the  respective  r  values  are  .45,  >50,  *40. 


The  ratio  of  Internal  energies  at  maximum  radius  and  Initial  radius 
is  E1/Eq  ■  (aQ/a1)^^"1^.  For  7  B  5/4,  these  ratios  are  0.077,  0.083# 
and  0.227  for  Pentolite,  RDX  +  Alum  (50/ 50),  and  2Hg  +  0g.  For  7  ■  4/3, 
the  ratios  decrease  to  O.O32,  0.035,  and  0.113. 

All  of  the  above  calculations  are  based  on  Friedman's  method.  If 
ve  use  a  least  squares  fit  of  the  experimental  data  to  Kennard's  eq.  I.30, 
ve  find  the  intercept  K  and  slope  KJa^  of  the  straight  line.  Using  the 
calculated  K  values,  and  assuming  1^  »  0.741,  the  respective  energy  par¬ 
tition  values  are  O.38,  0.4l,  and  0.33*  Using  the  calculated  KJa^  values, 
and  assuming  a^  =  1.0,  the  respective  energy  partition  values  are  0.49, 
O.67,  and  O.83.  (Assuming  a^  is  some  value  less  than  1.0  would  increase 
this  second  set.) 

If  Kennard's  eq.  1.30  were  a  better  fit  of  the  experimental  data 
than  Friedman's  eq.  1.23,  then  the  two  sets  of  data  should  agree.  Since 
they  differ  radically,  Friedman's  eq.  1.23  is  a  better  fit. 


CONCLUSIONS 


The  theoretical  equation  for  the  bubble  oscillation  period,  includ¬ 
ing  the  effects  of  surfaces,  was  found  to  be  Friedman's: 

t2  =  xvr^V5/6  (1.0  +  kh^V-^^p/Cd+b)) 

Using  the  method  of  least  squares  to  fit  the  experimental  periods 
of  three  selected  explosives,  the  period  constants  K  were  determined. 

The  error  in  K  was  calculated  to  be  less  them  5-5$  based  on  maximum  pro¬ 
bable  variations  in  the  experimental  variables.  The  energy  partition 
values  of  the  explosives  were  then  found  using  the  method  discussed  in 
Section  2.3*  This  involved  assuming  that  7,  the  ratio  of  specific  heats 
of  the  bubble  gas,  was  1.25.  The  resulting  r  values  for  Pentolite,  REX  + 
Alum,  and  2Hg  +  0g,  were  0.45,  0.50,  and  0.4l,  respectively. 

Increasing  7  would  decrease  the  resulting  energy  partition  values. 
For  example,  for  2Hg  +  0g,  if  7  is  increased  from  1.25  to  1.50,  r  is 
decreased  from  0.4l  to  O.38. 

Fran  Eq.  2.1,  the  accuracy  of  r  is  directly  proportional  to  the 
accuracy  of  Q  (probably  +  10$  at  best)  and  to  the  cube  of  K's  accuracy 
(+  17$) .  Thus  the  minimum  probable  experimental  accuracy  for  r  is 
+  27$,  and  the  energy  partition  values,  including  error  variation,  are 
0.45  +  0.12,  0.50  +  0.14,  and  0.4l  +  0.11,  respectively.  Although  the 
best  available  experimental  data  was  used,  the  minimum  probable  error 
variation  in  r  is  greater  than  the  variation  due  to  type  of  explosive. 

There  is  a  tendency  for  the  steam  bubble  to  retain  slightly  less 
energy  than  non-condensable  gas  bubbles  for  the  bubble  oscillation.  But, 
based  on  energy  partition  considerations,  it  makes  little  difference 
which  explosive  is  selected  to  simulate  the  nuclear  case. 
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APPENDIX  1 


THE  I  FUNCTIONS 

The  I  functions  are  used  in  determining  the  oscillation  period 
(1^  and  Ig),  the  vertical  momentum  (I^,  1^,  I^,  and  I g),  and  the 
vertical  displacement  of  the  bubble  (Iy  and  Ig).  They  are: 

I,  -  Z1  u-1  ^  da 

1  4o 

12  ■  }  u-1  a5'2  da 

13  -  £  u  a?/2  da 

IU  *  tv  *7/2  ^ 

15  ■  £  a'1  a’/2  ^ 

16  -  £  a’1  da 

I,  -  i  (a2  -  >2/2a h-1'*  a-3/2  da 

I6  -  £  (u2  -  3a2/2aV1/Z  a*1/2  da 

where  u  ■  (1  -  a^  -  J1/2. 

Graphs  of  the  first  six  I  functions  for  y  *  1*25,  1.33,  1.1*0,  l.!>0. 
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and  0  <  k  <  0.3,  are  given  in  reference  2.  A  check  on  their  .limit 
values  (at  k  ■  0),  using  the  beta  function 

1.0  , 

P(m,n)  ■  /  vra-  (l-v;11"1  dv  for  m  >  0,  n  >  0, 

0 

3 

where  v  -  a  ,  agrees  with  the  graphs  and  gives:  -  0. 71*68,  I2  -  0.6072, 

1^  -  0.1821,  -  0.1309,  Ig  -  0.1*668,  I6  -  0.1*250,  and  l*/  lg  -  0.9186. 
Except  for  Ij.  and  I^/Ig,  these  are  a11  upper  limit  values. 

This  method  cannot  be  used  to  determine  the  limit  values  of  1^  and 
Ig  since  they  are  integrated  from  a2  to  a^  and  a?  is  not  zero.  Evalu¬ 
ation  of  I-  and  Ig,  using  numerical  integration,  is  possible  but  very 
involved  because  of  the  s  term  (see  Appendix  2). 
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APPENDIX  2 


VERTICAL  MOMENTUM 


The  vertical  momentum  of  the  bubble  at  its  maximum  radius,  s^  ■ 
a^£^/3,  is  reported  in  reference  2  to  bet 

Bl  -  (L/Z)j^I5  +  LFI6/2(D+B)J  -  2L2F/(D+B)2jl3  -  LFI^/2(D+B)J 

where  F  *  dF/dx.  The  first  term  is  the  upward  momentum  due  to  gravity, 
and  the  second  term  is  the  downward  momentum  due  to  the  rigid  surface. 
The  vertical  momentum  at  the  bubble  minimum  is  twice  this  value. 

Usually  the  migration  at  maximum  radius  is  considered  negligible. 
However,  to  determine  the  effect  of  s^  in  Eq.  1.11,  the  calculations  in 
Table  A. 2  were  made. 

TABLE  A.  2 
Vertical  Momentum 


D  -  U.05  ft 

D  -  9.0  ft 

"l  ‘l/2  *1 

*1 

Pentolite 

-.0750 

HI 

♦.00002 

0.973 9 

RDX  ♦  Alum 

-.0696 

♦.002U 

0.9693 

+.0072 

♦.00003 

0.9717 

2H9  ♦  0o 

-.0200 

+.0002 

0.9200 

♦.0120 

+.00007 

0.9202 

Increasing  depth  increases  s^.  Since  s^  is  less  than  aero  at  shallow 
depth,  s12/2  decreases  to  a  negligible  value  at  mid-depth. 
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APPENDIX  3 
Table  A.  3 

Friedman's  Surface  Correction  Function:  F(x) 


X 

F(x) 

X 

F(x) 

0 

-  .05 

-  0.693 

-  0.785 

+  .05 

-  0.601 

-  .10 

-  O.878 

.10 

-  0.508 

-  .15 

-  0.975 

.15 

-  0.412 

-  .20 

-  1.076 

.20 

-  0.310 

1 

• 

TO 

VJl 

-  1.184 

.25 

-  0.202 

-  .30 

-  1.307 

•30 

-  O.O85 

-  .35 

-  1.431 

.35 

-  0.045 

-  .40 

-  1.577 

.4o 

0.191 

-  .45 

-  1.744 

.45 

0.358 

-  .50 

-  1.939 

•  50 

0.553 

-  .55 

-  2.174 

.55 

0.788 

-  .60 

-  2.462 

.60 

1.076 

-  .65 

-  2.829 

.65 

1.443 

-  .70 

-  3.312 

.70 

1.926 

-  .75 

-  3.985 

.75 

2.599 

-  .80 

-  4.991 

.80 

3.&5 

-  .85 

-  6.661 

.85 

5.275 

-  .90 

-  9.998 

.90 

8.612 

-  .95 

-20.000 

.95 

18.620 

-1.00 

-  00 

1.00 

00 

The  following  series  expression  was  found  to  agree  with  tbs 
above  tabulation: 

f  -  -0.693  +  1.832x  +  1.978X3  +  ac5(l-x2)-1 
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APPENDIX  U 


EXPERIMENTAL  DATA 


Three  explosives  were  studied  in  the  Hydra  test  pond  during  the 
summer  of  I960.  They  were:  Pentolite,  the  control}  an  equal  mixture 
by  weight  of  HDX  +  Alum  (A1  NH^(S0^)2*12H20)}  and  a  pressurized  stoichi¬ 
ometric  mixture  of  hydrogen  and  oxygen  in  a  frangible  plastic  sphere. 

All  charges  were  spherical  and  centrally  detonated.  They  were  detonated 
at  various  depths  in  the  hemispherical  pond  of  16-ft  radius.  Bubble 
radii  and  periods  were  measured  • 

MAXIMUM  BUBBLE  RADIUS 

A  high  speed  Photosonic  camera  was  mounted  underwater,  about  25  ft 
from  the  charge  (in  a  camera  bay).  A  six  inch  wire  grid  was  mounted 
midway  between  the  camera  and  the  charge.  Films  showing  a  well  defined 
first  oscillation  were  reduced  using  the  grid  lines,  for  the  first  maxi¬ 
mum  horizontal  diameter.  The  experimental  J  values  “  az^At1/3) 

are  compared  with  the  theoretical  J'a^  values  in  Table  A.lwl. 


TABLE  A.U.1 
Radius  Constants 


Explosive 

j,ai 

Pentolite 

12.96  i  0.1*8 

13.07 

RDX  +  Alum 

13.01  *  0.2U 

13.U7 

2H2  ♦  02 

11.99  *  0.5U 

11.97 

*  Complete  experimental  details  can  be  obtained  from  W.W.  Perkins 
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BUBBLE  PERIOD 

Bubble  periods  were  measured  by  recording  on  an  oscilliscope  camera, 
the  pressure  signal  from  a  tourmaline  piexoelectric  gauge  3et  about  12  ft 
from  the  charge.  Initially  the  scope  swept  twenty  milliseconds  per  centi¬ 
meter,  recording  the  time  between  the  shock  wave  and  the  bubble  pulse. 
Later  in  the  series,  in  order  to  more  accurately  determine  the  depth  ef¬ 
fect  on  the  period,  the  scope  was  delayed  a  set  time  (150  to  180  msec) 
from  the  initiating  pulse  and  then  swept  slower  (from  0.5  to  5«0  msec/ 
cm)  to  record  the  bubble  pulse  only.  In  some  cases  the  estimated  delay 
time  was  incorrect  and  the  bubble  pulse  was  partially  or  entirely  missed. 
But  where  the  period  was  obtained,  it  was  at  least  four  times  more  accu¬ 
rate  than  those  from  the  earlier  method. 

The  following  tabulated  data  are  periods  obtained  using  the  delay 
method.  Since  the  same  pulse  triggered  both  the  explosive  and  the 
scope  delay  circuit,  the  observed  time  on  the  oscilliscope  included  the 
time  of  arrival  of  the  first  minimum  pressure  pulse.  This  arrival  time 
was  found  to  be  2. It  msec  for  Pentolite  and  BOX  +  Alum  (50/50),  and  2.8 
msec  for  2Hg  +  Og*  Tha  Periods  1x1  Tables  A.U.2  -  A.h.i*  have  been  cor¬ 
rected. 
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table;  a.4.2 


Experimental  Pentolite  Periods 
WEIGHT:  1.06  lb 


GAUGE:  0.5  inch  diameter,  except  vhere  noted 

DETONATOR:  Engineer's  Special  (0.875  g  PETN,  plus  mercury  fulminate) 


Shot 

D 

(ft) 

D+B 

(ft) 

T 

exp 

(sec) 

T 

calc 

(sec) 

AT 

(sec) 

66 

4-3 

18 

.1906 

.1850 

-.0056 

61 

4.3 

18 

.1896 

.1850 

-.0046 

67 

4.3 

18 

.1896 

.1850 

-.oo46 

50 

5.0 

18 

.1866 

.1892 

+.0026 

51 

5.0 

18 

.1875 

.1892 

+.0017 

901 

5.0 

18 

.1859 

.1892 

+.0033 

901* 

5.0 

18 

.1854 

.1892 

+.0038 

902 

5.0 

18 

.1851 

.1892 

+.oo4i 

902* 

5.0 

18 

.1842 

.1892 

+.0050 

75 

7.0 

19 

.1956 

.1936 

-.0020 

74 

9.5 

19 

.1966 

•  1937 

-.0029 

*1.0  inch  diameter  gauge. 
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TABLE  A.4.3 

Experimental  REK  +  Alum  (50/50)  Periods 

WEIGHT:  0.8  lb  of  RDX,  plus  0.8  lb  of  Alum 

GAUGE:  0.5  inch  diameter,  except  vhere  noted 

DETONATOR:  Modified  Engineer's  Special  (0.875  grams  tetryl,  plus 
lead  azide,  plus  a  2.5  gram  tetryl  booster) 


Shot 

D 

(ft) 

D+B 

(ft) 

T 

exp 

(sec) 

T 

calc 

(sec) 

AT 

(sec) 

90* 

4.0 

18 

.1746 

.1749 

+.0003 

96* 

4.0 

18 

.1770 

.1749 

-.0021 

96 

4.0 

18 

.1772 

.1749 

-.0023 

94* 

5-0 

18 

.1792 

.1818 

+.0026 

94 

5.0 

18 

.l8o4 

.1818 

+.0016 

93* 

6.0 

18 

.1826 

.1854 

+.0028 

56 

9.5 

19 

.1882 

.1871 

-.0011 

57 

9-5 

19 

.1886 

.1871 

-.0015 

*1.0  inch  diameter  gauge. 


TABUS  A.4.4 

Experimental  2Hg  +  Og  Periods 

GAUGE:  0.5  inch  diameter,  except  where  noted 
DETONATOR:  Pyrofuze  (Ra.llad.ium  and  Aluminum  alloy  wire) 


Shot 

w 

(H>) 

D 

(ft) 

D+B 

(ft) 

T 

exp 

(sec) 

Tcalc 

(sec) 

AT 

(sec) 

62 

0.3175 

4.0 

19 

.1712 

.1&9 

I 

• 

8 

S' 

70 

0.2635 

4.0 

18 

.1572 

.1581 

+.0009 

77 

0.2635 

6.0 

19 

.1582 

.1683 

+.0101 

87* 

O.2635 

6.0 

18 

.1672 

.1686 

+.0014 

92* 

O.2635 

6.0 

18 

.1078 

.1686 

+.0008 

92 

O.2635 

6.0 

18 

.1687 

.1686 

-.0001 

l6 

0.2635 

8.0 

19 

.1732 

.1707 

-.0025 

6l 

0.3175 

9.5 

19 

.1862 

.1812 

-.0050 

*1.0  Inch  diameter  gauge. 
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APPENDIX  5 


METHOD  OF  LEAST  SQUARES 


The  best  linear  fit  to  a  series  of  values  is  a  line  about  which  the 

sum  of  the  squares  of  the  deviations  is  a  minimum.  Applying  this  prin- 
2 

cipal  to  y  -  K  +  HK  z,  and  minimizing  the  error  in  y,  we  findt 

K  -  ($yEz*  -  Zz£yz)/(N£z2  -  EzZz) 

HK2  ■  (EyEz  -  N?yz)/(£zEz  -  NEz2) 

where  N  ie  the  number  of  data  sets  (y,z). 

If  the  error  in  z  were  minimized,  then: 

K  ■  (Ey?yz  -  £z$y2)/(N2yz  -  Zylz) 

HK2  -  (Sy£y  -  N2y2)/(HjrEz  -  Ntyz) 
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APPENDIX  6 


ERRORS  IN  K  AND  H 


From  the  method  of  least  squares,  K  is  proportional  to  y  which  equals 
and  H  is  proportional  to  l/yz  which  equals  (D+Bj/TFZ^^.  Thus 
any  errors  in  the  measurement  of  the  period  and  or  the  depths  will  affect 
the  accuracy  of  K  and  H.  K  will  also  be  affected  by  weight  errors.  From 
the  experimental  data,  the  period  errors  and  the  weight  errors  are  each 
less  than  2%.  To  determine  the  effect  of  the  depth  error,  the  following 
assumptions  were  made:  (l)  the  error  in  charge  depth  D  due  to  water  sur¬ 
face  wave  motion  is  less  than  0.18  feet,  and  (2)  the  pond  depth  error, 
due  to  water  blown  out  by  preceding  shots,  is  less  than  0.5  feet. 

Then  for  this  data,  the  0.5  foot  variation  in  pond  depth  will  cause 
less  than  U»7%  error  in  (D+B)/F.  The  0.18  foot  variation  in  charge  depth 
will  cause  less  than  0.1$  error  in  Z^  and  less  than  5.6$  error  in  FZ1/2. 
The  hydrostatic  charge  depth  in  sea  water  is  D+33  feet  while  in  fresh 
water  it  is  D+3U  feet.  The  variation  is  less  than  2.356  in  Z^6  and  less 
than  1.1$  in  Z^2. 

Thus,  from  the  above  assumptions,  the  total  error  in  K  is  less  than 
5.5%  and  the  total  error  in  H  is  less  than  11$. 

An  additional  possible  error  in  the  experimentally  determined  K 

•it 

should  be  mentioned  here.  Rudlin's  work  indicates  that  for  spherioal 
explosives,  the  fraction  of  the  charge  consumed  is  a  function  of  the 


*  L.  Rudlin.  An  Approximate  Solution  of  the  Flow  Within  the  Reaction 
Zone  Behind  a  Spherical  Detonation  Wave  in  TNT,  U.S.  Naval  Ordnance 
Laboratory.  NAVWEPS  7361;,  April  1961. 
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charge  size.  Assuming  Bullin' s  theory  for  TNT  applies  for  Pentolite, 
the  following  values  were  calculated.  Here  N  is  the  fraction  of 
explosive  consumed. 


W(lbs) 

N 

0.05 

0.759 

1.06 

0.880 

106.00 

0.942 

The  experimental  value  of  K1  was  found  to  be  4.30  for  a  1.06  pound 
Pentolite  charge.  But  if  only  88£  of  the  explosive  was  consumed,  the 
actual  weight  was  .88  x  1.06  pounds.  Since  is  constant,  the  ex¬ 

perimentally  determined  K'  value  can  be  predicted  to  vary  with  the 
charge  weights  K'  ■  4*49  ' .  In  other  words,  the  K'  value  found 

wp 

with  a  1.06  pound  Pentolite  charge  should  be  2%  less  than  the  K1  value 
found  using  a  416  pound  charge  and  5%  higher  than  the  K1  value  found 
using  a  0.05  pound  charge.  This  effect  should  not  apply  to  the  SHg+Og 
explosive. 
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APPENDIX  7 

3YMB0L3 


All  dimensions  are  in  the  fps  system 
A  bubble  radius 

a  non-dimensional  radius 

B  charge  distance  from  bottom 

b  B/L 

C  unit  of  time 

D  charge  depth 

d  D/L 


E 

F,  F(x) 
f(x) 

g 

H 

I 

J 

K 

KE 

k 


internal  bubble  energy 

function  giving  surface  effects  on  period 

function  giving  surface  effects  on  period 

acceleration  of  gravity 

second  order  bubble  period  constant 

see  Appendix  I 

bubble  radius  constant 

bubble  period  constant 

kinetic  energy  of  water 

a  3(r~ i) 
o 


29 


L  unit  of  length 

P  bubble  pressure 

Pk  hydrostatic  pressure  at  charge  depth}  ■ fgZ 

Q  charge  energy  per  unit  weight 

r  fraction  of  charge  energy  in  bubble 

s  non-dimensional  vertical  momentum  of  bubble 

T  time  variable 

Tg  time  to  first  raimimum,  i.e.  bubble  period 

t  non-dimen3ional  time  variable}  »  T/C 

V  bubble  volume 

W  charge  weight 

WD  work  of  displacement 

x  (0-B)/(D+B) 

Z  D+33  ft  (3U  for  fresh  water) 

Y  ratio  of  specific  heats  of  bubble  gas 

p  mass  density  of  water 
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